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Abstract 

Within the Euclidean path integral and mass perturbation theory we derive, from the Dyson- 
Schwinger equations of the massive Schwinger model, a general formula that incorporates, for 
sufficiently small fermion mass, all the bound-state mass poles of the massive Schwinger model. 
As an illustration we perturbatively compute the masses of the three lowest bound states. 
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1 Introduction 

The massless Schwinger model - which is two-dimensional QED with one massless fermion - is 
wellknown to be exactly soluble (PJ - {§, |H] - pO(j), and its solution may be used as a starting 
point for a (fermion) mass perturbation theory of the massive Schwinger model (]l^|, 118(1 , iP^ li 
]l9|) . In both models instantons and a nontrivial vacuum structure (^-vacuum) are present (0 
- fl2|] ). The spectrum of the massless model consists of one free, massive boson with Schwinger 
mass /Up = (fermion-antifermion bound state ||, ||) and trivial higher states consisting of n 
free Schwinger bosons. 

For the massive model these higher states turn into n-boson bound states. Their masses, 
in principle, could be computed, using mass perturbation theory, by evaluating the mass poles 
of the corresponding n-point functions. Here we will adapt a slightly different method. By 
exploiting the Dyson- Schwinger equations of the model we will find that all bound-state mass 
poles are contained within one formula. From this we will compute the masses of the three 
lowest bound states perturbatively. 

All computations are performed within the Euclidean path integral formalism and are done 
for general vacuum angle 0. This latter fact causes some minor complications, because for 8^0 
parity is no longer conserved and, as a consequence, the mass pole equations will turn into 
matrix equations. 



2 Massless Schwinger model 

Before starting the actual computations, we need some formulae from the massless Schwinger 
model. Indeed, the vacuum functional and Green functions of the massive Schwinger model 
in mass perturbation theory may be traced back to space-time integrations of VEVs of the 
massless model, 

oo 

Z(m,6)= e ik9 Z k (m) (1) 

k=— oo 

{k . . . instanton number) where 

„ oo n n /• 

Z k {m) = N J D^D^DA" k £ — J[ J ^(x^^)- 



/ dx 

e J 



(2) 



and the mass perturbation expansion is yet performed. A general VEV of the massive model is 
given by 

{0)m = Z(^9) {0 £ IT / dxM^M^W (3) 

^ ' ' n=0 ' i=l 

For these expressions we need (pseudo-) scalar and vectorial VEVs of the massless model. It is 
useful to rewrite the scalar densities in terms of chiral ones, S(x) = S+(x)+S-(x), S± = *f>P± i f?, 
because for VEVs of chiral densities only a definite instanton sector k = n + — n_ contributes, 
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(S Hl ( Xl ) ■ ■ ■ S Hn (x n )) = e ike (^Y exp&i-r^Dnixi - x 3 ) 



(4) 



i<3 



(see e.g. Q^Hi [pOj , Jjj]] for the computation) where a% = ±1 for H L = ±, D^ is the massive 
scalar propagator, 

d im>( x ) = -7T K o(^o\x\), Dpoip) = -s- — 2' ( 5 ) 



(i^o • • • McDonald function) and S is the fermion condensate of the massless Schwinger model, 



S = (M) = — ju 

27T 



(6) 



(7 . . . Euler constant). An inclusion of an arbitrary number of vector currents does not alter the 
contributing instanton sector and may be computed from the generating functional 

(S Hl ( Xl ) ■ ■ ■ S Hn (x n M(3] = e^(|)"exp[^(-)^^4^ (x i - Xj ) 



Kj 



■ exp 



dyidy 2 (5{y 1 )D^{y 1 - y 2 )f3(y 2 ) + 20F^(-)^ / dyf5{y)D^{y - x t ) 



1=1 



(7) 



More precisely, (7) generates all VEVs of n chiral densities and an arbitrary number of Schwinger 
bosons 6, where 6 is related to the vector current via 



(8) 



(7) may be found by the inclusion of a vector current source into the path integral quantization 
and was explicitly computed in pi]. 



3 The bound-state mass poles 

First we have to fix some notation for later convenience: 

E±(x) := e ±4 ^W*) _ l 



1=0 



E^\p) = J d 2 xe^ x E^\x) , 4 n) :=E^ ] (0). 



(9) 



We will use the following Feynman rules: 
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m 

go 
1 



vAAAA/V 



X 



Fig. 1 



where m and go = m(S(x)) m are the bare and renormalized coupling, respectively. In the 
sequel all VEVs are with respect to the massive model, therefore we will omit the subscript m. 

We will discuss the special case 9 = first, because it is easier and may be represented by 
simple graphical computations. Later we generalize to arbitrary 6. 

On fermionic bilinears there hold two equations of motion, namely the Maxwell equation 
and the anomaly equation. Eliminating the field strength one arrives at 



(□ x - f4)<t>(x) = 2vWP0r) 



P 



(10) 



where the Schwinger boson (f> is related to the vector current like in (8). Introducing the abbre- 
viation 



(11) 



one may derive Dyson- Schwinger equations like e.g. for the two-point function, 

My.My^yMyi)) = M yi S( yi - y 2 )+ 
47T 50 <%1 - V2) + 47T^(P(yi)P(y 2 )), 



(12) 



where it is understood that external sources couple with coupling constant 1. The validity of 
equation (12) is most easily seen in a graphical representation for the two-point function, 



+ 



+ 



+ 



+ 



Fig. 2 
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Here all graphs where the two boson lines meet on one point contribute to the renormalized 
coupling, 



+ S + S + ... 



Fig. 3 



g = m(S) = mS + ^(m£) 2 (P+ +£_) + ..., (13) 



see |L7j], [18] for computational details. All internal vertices are renormalized by the same cou- 



pling g , and the remaining (P(yi)P(y 2 )) part may be written like 



+ 



Fig. 4 

The {P{y\)P{y 2 )) propagator includes, even in least order, an arbitrary number of bosons 
propagating from y\ to y 2 , which will be essential in the following. 

Similar Dyson-Schwinger equations may be derived for higher n-point functions, e.g. for the 
four-point function 

M yi M y2 My 3 M yA {4>{ yi )(j){y 2 )4>{y^{y^) = 

M yi M y2 {cj)(y 1 )(P(y 2 ))M y3 M yi {<j)(y 3 )(P(y 4 ))+ perm. + 

I6ir 2 g 5(yi - y 2 )5(yi - ysjSfa - y 4 )+ 

16vr 2 5fo5(yi - 2/2)<%3 - y4){S(yi)S(y s )) + perm. + 

16vr 2 9 ^(yi - y 2 )(S(y 1 )P(y 3 )P(y i )) + perm. + 

167r 2 g^P( yi )P(y 2 )P(y 3 )P(y A )) (14) 

and analogously for higher n-point functions. 

The essential point is that in all these Dyson-Schwinger equations there occurs an identical 
term that will be responsible for the bound-state formation for sufficiently small fermion mass. 
In momentum space this term reads 

c(g +g 2 (PP)(p)) (15) 
for odd bound states (c = 4-7T, 167r 2 , . . .) and 

c(g + g 2 (SS)(p)) (16) 
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for even bound states. 

Now both terms (15), (16) may be inverted via the geometric series formula, e.g. 

go(l+go(PP)(p)) = y ~ - (17) 

l-g (PP) ni (p) 

where n.f. stands for nonfactorizable and means that graphs contributing to (PP) n f may not 
be factorized in momentum space. Its graphical representation looks like 




XwX + 



Fig. 5 

In (17) the small <?o has to be compensated by a large contribution from (PP) n{ in order 
to give rise to a mass pole. The first term in Fig. 5 is E±(x) = e ±47rDA1 o ( x ) — 1 and contains 
an arbitrary number of massive propagators D^ a . Now precisely (D fJi0 (x)) n has a threshold 
singularity in momentum space at p 2 = (rifio) 2 , therefore (17) may have mass poles near p 2 = 
(n^ ) 2 for n £ N. More precisely, rewriting (for = 0) {PP) = 2(S + S + ) - 2(S + S_) (and with 
a + for (SS)), and using 

(S + (x)S + (0)) = (|)V 4 ^oM j ( S+ (x)S-(0)) = (|) 2 e - 4 ^oW (is) 

we find that (PP) may cause mass poles for odd n whereas (SS) may cause mass poles for even 
n, as it has to be. For all mass poles only the terms (15), (16) may balance the pole equation 
for sufficiently small go> therefore it is enough tojxtnsider them. 

To get more insight we next have to rewrite (PP) n f (Fig. 5) in terms of internal Schwinger 
bosons (we ignore constants) 

+ ... 



Fig. 6 

We find that the one-boson propagator acquires no corrections, whereas all terms with more 
than one boson propagating from y\ to yi have corrections along the boson lines. Consequently, 
one can compute the lowest pole mass (the Schwinger boson mass) without the need to know 
it. 
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On the other hand, for the computation of higher bound-state masses, one has to take into 
account the corrections, i.e. insert the exact Schwinger mass (the lowest pole mass). The reason 
is that the mass corrections for the bosons just shift the position of the threshold singularity 
and are therefore important in lowest order. There are other corrections present, too (internal 
boson interactions), however, they are unimportant in lowest order. 

This result is very plausible physically: the higher bound states should consist of physical 
Schwinger bosons with their physical masses \x (not the bare masses n§). 

All these features remain true for general 6, only the pole mass formula itself is slightly more 
complicated and shall be derived next. 

4 General 9 case 

For 9 7^ the renormalized coupling is complex, 

g = m(S + ) , g* = m(S-) (19) 

and the Dyson-Schwinger equations are slightly changed, too. E.g. the two-point function obeys 

M yi M y2 {4>{ yi )cj){y 2 )) = M yi S{ yi - y 2 ) + (g + g*)6( yi - y 2 )+ 

g 2 (S + ( yi )S + (y 2 )) + (g*) 2 (S-( yi )S-( y2 )) - 2gg*(S + {y 1 )S-(y 2 )). (20) 
The interesting function that gives rise to the mass poles is 

g(l+g(S^S + )( P ) - g*(s7S-)(p))+g*(l + g*(sZS-)(p) - g(S^S-)(p)) (21) 

for odd bound states and with only plus signs for even bound states. Because parity is no longer 
conserved the P and S components mix and the geometric series formula (17) generalizes to a 
matrix equation. 

Introducing the abbreviations 

A:=l+g(s7s + )(p)-g*(s7s^(p) 
a:=g(sTs + ) Q , L (p) , P ■= g(S^S-)^(p) (22) 

the equation reads 

A = 1 + aA- (3* A* 

A* = 1 + a* A* - 0A (23) 

and has the solution 

A = — . (24) 

1 - a - a* + aa* - [3(3* K J 

Equation (23) may be checked by a careful investigation of the perturbative expansion for 
(P( yi )P(y 2 )). 

For even bound states the solution may be found from (24) by the substitution f3 — > — (3. 
Therefore, all solutions have the same denominator, and the zeros of this denominator are the 
pole masses of all bound states of the massive Schwinger model. 
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Explicitly the pole-mass equation reads 

(1 - g(S^S+) luL ( P ))(l - g*(s7s + )* n Ap)) = 99*((s7S-) n . L (p)) 2 (25) 

where 

g = m ^e ie + m 2 (^) 2 (E + e 2ie + E_)+o(m 3 ) 

= ■ 9i+92+o(m 3 ) (26) 

(sTs + ) n Ap) = E+(p) + o(m) , (sTS-) n , L (p) = E-(p) + o(m). (27) 
Of course, for 9 = (real A,a,{3), one recovers the equations (15), (16). 

5 Explicit mass computations 

For a computation of the Schwinger mass up to second order we rewrite (25) like 

(1 - gE+(p))(l - g*E + (p)) = gg*E 2 _(p) (28) 
and separate the one-boson contribution 

E±(p) = -4^ + E±(p) (29) 

leading to 

- p 2 - i4 = SirRe g x - 2Re g x Ef (p) (p 2 + p 2 ) + SvrRe g 2 + 8ir 9l gl (E^ (p) - E W (p)) (30) 
with the solution 

-p 2 = p 2 1 (l + 4 7 r 



p 2 = nl(l + 47r-^cos6<) (31) 



in first order and 

Em „ m 2 E 2 



2 2 

P = t*o 



1 + 4^— cos + 2^ — — +^ 1) (l))cos26' + J E;_ -^(l))] (32) 
p p 



in second order. Here we rescaled p — > p' = and used EyP (p') = E± (l)+o(m) in the last step. 
This result precisely conicides with the result obtained by a direct perturbative computation 
(0)- 

In order to compute the two-boson bound state we have to separate the two-boson part of 
E+ in (28). In lowest order we find 

1 = ^(5i + ^)16tt 2 (^)(p) (33) 
where now p is the physical Schwinger mass (32) including fermion mass corrections. Using 
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(see e.g. [22[ for a computation) and remembering that 4p 2 — (— p 2 ) is a very small number 
(that is positive for a bound state) we may set — % — — arctan(- • •) ~ ^ and get 



_4 rrl 2y2 

p 2 ^V(l-^-^cos^) 
16 p 4 



(35) 



which is of second order in m. Again, this result coincides with the one from a direct perturbative 
calculation (l2l 



6 The three-boson bound-state mass 

For the three-boson bound-state mass we have to separate the three-boson part in (28) and 
find, in lowest order 

1 = ^mEcosfl • 64vr 3 (L>3)(p) (36) 
or, after a rescaling p — > ^ to dimensionless momenta 



64vr 3 mS ~ 
1 = — -cos0(Z*>)(p). 



(37) 



6 /i 2 

is given by the graph (where we introduce positive squared momentum Q 2 = —p 2 > 0) 




Fig. 7 



d 2 q 1 d 2 q 2 



1 1 



(2vr)4 (p + gi + (6 )2 + lg2 + lg2 + 1 



where 



-2 1 dx I dy I 
Jo Jo J 



d 2 q 1 d 2 q 2 



1 



(2tt)< 



dx 



Qi + 1 + (?2 - ?i> + ((p + 9i + ^) 2 - g!)j/ 



(47r) 2 7o Q 2 (xy — x 2 y — y 2 + xy 2 ) — x + x 2 — xy + y 2 
g?x ft d-Z 



/ 



8vr 2 (l-Q 2 (l-x))7 z 2 + T 2 (Q 2 ,x) 



fix 



1 X 

■ arctan 



o 8tt 2 (1-Q 2 (1-x))T(Q 2 ,x) 2T(Q 2 ,x) 
x 2 - Q 2 x 2 (l - x) + 4x(l - x) 



T 2 (Q 2 ,x) 



4(Q 2 (l-x)-l) 



(38) 
(39) 
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12vr2 7 y/\9x-8\ /( x _ *)2 X + ex 2 (1 



The numerator of T 2 has a double zero at Q 2 = 9: 

9x(x - ^) 2 . (40) 

This double zero is in the integration range of x and is precisely the threshold singularity. 
Setting 

Q 2 =: 9(1 - e) (41) 

in the numerator of T 2 in the factor A, and Q 2 = 9 everywhere else, where it is safe, one arrives 
at: 

^/x\9x-8\ 

j arctan — 7— 

<h 3{X " 3) ==:I(e). (42) 

+ ex' z (l — x) 

The mass-pole equation reads (^cos# = ^2-2Re<?i =: a) 

647T 3 

1 = ^-al(e) (43) 

and must be evaluated numerically. It gives rise to an extremely tiny mass correction e. For 
sufficiently small a it is very well saturated by 

e(a) - 0.777 exp(-^^) (44) 
a 

and is therefore smaller than polynomial in a. (I checked the numerical formula (44) for 30 < 
^ < 1000, corresponding to 10~ 3 < e < 10 -100 , but I am convinced that it remains true for 
even larger — ; however, there the numerical integration is quite difficult because of the pole in 
(42).) 

A more accurate mass-pole equation would include some additional contributions: 

647T 3 

l = a(— 1(c) + /(c)) (45) 

where /(e) is some function that is finite for e — ► 0. But for a sufficiently small it remains true 
that an extremely tiny value of e suffices to saturate the mass-pole equation, whatever the value 
of /(e) is. 

We conclude that the three-boson bound state mass is nearly entirely given by three times 
the Schwinger boson mass, 

Mf ~ 9/x 2 , (46) 

or, differently stated, that the binding of three bosons is extremely weak. 

This result enables us to add a short remark on a result that was obtained in There it 
was argued that, for 9 = 0, the three-boson bound state should be stable because a decay into 
two Schwinger bosons is forbidden by parity conservation. However, because of (46) it holds 
that M3 > M2 + (jl, therefore there should be a small probability for the three-boson bound 
state to decay into one two-boson bound state and one Schwinger boson. 
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7 Summary 

As claimed, we arrived at our aim to derive one formula for all bound-state masses of the 
massive Schwinger model, at least for sufficiently small fermion mass. Of course, if we were 
able to exactly solve the model, it would not at all be surprizing to find all mass poles within 
one Green function. The interesting point is that we could reach this aim by the use of the 
Dyson- Schwinger equations and by a specific partial resummation of the perturbation series. 

We found that the two lowest states acquire noticeable corrections, whereas the binding 
energy of the three-boson bound state is extremely tiny. It is plausible to conjecture that the 
binding energies of higher bound states remain very tiny. 
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